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Abstract. We consider DNA-cationic lipid complexes that form lamellar stacks of lipid bilayers with parallel 
DNA strands intercalated in between. We calculate the electrostatically induced elastic deformations of the 
lipid bilayers. It is found that the membranes undulate with a periodicity that is set by the DNA interaxial 
distance. As a consequence the lamellar repeat distance changes resulting in a swelling or compression of 
the lamellar stack. Such undulations may be responsible for the intermembrane coupling between DNA 
strands in different layers as it is observed experimentally. 

PACS. 68.10.-m Fluid surfaces and fluid-fluid interfaces - 64.70.Md Transitions in liquid crystals 



1 introduction 



Electrostatic adsorption of polyelectrolytes onto oppositely 
charged surfaces, such as lipid membranes, has been the 
subject of intense experimental and theoretical research 
in the last decade. Of particular interest is the sponta- 
neous complexation of DNA with both cationic and neu- 
tral lipids due to their possible application to gene ther- 
apy I. These so-called "lipoplexes" show a diversity 
rium and metastable structures ^,^,||,||,|^,||,||, 
12 1 . For example, it has been shown through X-ray 
0,0 that DNA molecules and hpids 



of equi 



diffraction analysis 

can form lamellar complexes with DNA intercalated in 
between lipid bilayers. Another complex formed for suf- 
ficiently flexible membranes is an inverse hexagonal lipid 
structure with DNA inside the water regions m. 

Several theoretical studies help to understand many of 
the phenomena observed for lipoplexes |]l3, 14, r^,[l6[^, 
|l|,|l9|,|20[|2^,||,||,|24[|2|,|^,|3, however many more re- 
main to be elucidated. May et al. studied in detail the 
phase behavior of aqueous mixtures of DNA, cationic lipid 
and neutral lipid. Their model is based on the two ther- 



which the anionic charges of the DNA balance the cationic 
charges of the lipids) is unstable to further adsorption of 
DNA or lipids. The formation of say an isoelectric com- 
plex is driven by the release of the small counterions that 
were " condensed" on the highly charged DNA and on the 
charged bilayer before complexation. A lipoplex close to 
the isoelectric point is very susceptible to the uptake of 
further cationic lipids or DNA ~ if available - since this 
will be accompanied by the release of the corresponding 
counterions into the lipoplex]^. The theoretical predictions 
show good agreement with a recent experimental study 

A different approach to the problem of the distance be- 
tween DNA strands was given by Dan [Q. In this study 
the preferred distance was predicted to be the result of two 
competitve mechanisms, electrostatic repulsion between 
the strands and their membrane-induced attraction due 
to the perturbation of the lipid packing in the membrane 
close to the adsorbed DNA. The model assumes DNA 
strands adsorbed on a single membrane as it was inves- 
tigated experimentally by Fang and Yang Q using atomic 
force microscopy experiments. They found that the dis- 
tance between DNA strands was about 5 nm, a distance 



modynamically stable structures found experimentally fcii , i , j • j. j u t-v -j-u ■ j i ii 

,. r ;• r;i T^T^TAn- . that was predictcd by Dau withm hcr modcl 1 

and predicts the phase diagram as a function oi the DN A/lipid r ^ l 

composition and the elastic properties of the lipid bilay- 
ers. Another problem that was considered theoretically is 
the dependence of the interaxial spacing of DNA rods in 
lamellar complexes on the DNA/lipid composition. Bru- 
insma 1 18 presented an analytical approach that is appli- 
cable to lipoplexes with weakly charged bilayers. The nu- 
merical study of Harries et al. predicts the interaxial 
spacing also for higher charge densities. According to both 
studies the isoelectric point of the lipoplex (the point at 



A similar instability is also expected for the complexation 
of charged spheres and oppositely charged polyelectrolytes. A 
single highly charged chain will wrap around a single sphere 
forming a complex that is beyond the isoelectric point ("over- 
charging" ) and this effect is driven by the release of counterions 
from the wrapped chain j2^. On the other hand, a chain in a 
solution of highly charged spheres will complex more spheres 
than necessary to be isoelectric, and this effect is driven by the 
release of counterions of the complexed spheres [E9|. 
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In that model it is assumed that the membrane (that 
is supported by a soHd surface) is locally perturbed close 
to the DNA in such a way that the monolayer thickness 
is slightly increased jl^. It should be expected that for 
lamellar lipoplexes one also has perturbations. Since the 
membranes are allowed to undergo shape changes freely 
(no supporting layer) one might expect undulations lead- 
ing to a compression or sweeling of the whole lamellar 
stack as depicted in Fig. 2. Such undulations might lead 
to an intermembrane coupling between DNA rods in dif- 
ferent layers - resulting in a 3D ordering of the DNA rods. 
Lipoplexes with a 3D rectangular ordering of the DNA 
molecules were indeed observed experimentally Q . 

Most experiments on lamellar lipoplexes indicate that 
such a type of perturbation of the membranes around the 
DNA molecules - if present - is small [|[|^. Undulations of 
the membranes should lead to a lamellar repeat distance 
that is larger or smaller than the sum of the bilayer thick- 
ness and the diameter of the DNA molecule (including a 
hydration shell). Considering complexes at the isoelectric 
point and changing the ratio of charged to neutral lipids 
it was observed that the lamellar repeat distance stays 
always close to a value that indicates flat membranes g 
1^. Thus even though the lipid dilution experiments lead 
to a considerable increase of the interaxial spacing be- 
tween DNA rods, the undulations remain too small to be 
non-ambiguously detected. On the other hand, for more 
flexible membranes where detectable membrane undula- 
tions could be expected the system switches to the inverse 
hexagonal phase instead . 

Recently Subramanian et al. [ ^0[ studied the complex- 
ation of the anionic polypeptide poly-glutamic acid with 
a mixture of cationic (DDAB) and neutral (DLPC) lipids 
by means of small angle X-ray scattering and neutron 
scattering. It was observed that the lipid organizes in a 
multilamellar phase with the polypeptide chains interca- 
lated in between the membranes. Compared to the DNA 
complexes discussed above, the polypeptides do not show 
any in-plane ordering even though it is assumed that they 
are in the a-helical state. As for the DNA lipoplexes a 
"lipid dilution" experiment was performed for isoelectric 
polypeptide lipoplexes. Contrary to the outcome for the 
DNA complexes, a considerable increase of the lamellar 
spacing was found when the cationic lipids were diluted by 
neutral ones. For high lipid dilution the spacing saturated 
at a constant value of 60A which coincides with the equi- 
librium value of pure DLPC membranes. Subramanian et 
al. suggested that this behavior could be due to a 
"pinching mechanism" including membrane undulations 
similar to the ones depicted in Fig. 2 (case h < 0). The 
pinching sites are formed due to the electrostatic inter- 
action between the negatively charged poly-glutamic acid 
and the cationic DDAB lipids. Away from the pinched re- 
gions the properties of the lipoplex are dominated by the 
properties of the pure DLPC membranes. Whether it is 
possible to have pinches in a lipoplex was studied by one 
of the authors p2| . By comparing the gain in electrostatic 
free energy with the bending energy of forming a pinch, 
the parameter range was estimated at which pinching can 



be expected. It was shown that this effect should occur if 
the line charge density of the rods is sufficiently high and 
the membranes are sufficiently flexible, a situation that 
might be fulfilled for the polypeptide lipoplex considered 
in Ref. @. 

A different approach to the pinching problem is taken 
in the present study. We start out with a perfectly flat 
lamellar lipoplex as depicted in Fig. 1. The DNA rods are 
assumed to be ordered within a 3D rectangular lattice as 
it was observed by Artzner et al. Q. Our goal is to cal- 
culate how the electrostatic interaction between the neg- 
atively charged "rods" and the positively charged mem- 
branes modifies the conformation of the membranes. We 
show that there are in principle two possibilities, namely 
a compression of the lamellar stack as depicted in Figure 
2 (/i > 0) or an expansion as depicted in the same Figure 
(case h < 0). 

In the next section we introduce the model system and 
calculate its electrostatic and bending free energies for ar- 
bitrary but small periodic undulations of the membranes. 
By minimizing the free energies of the undulation with 
respect to its Fourier components we show in Sect, ^that 
the electrostatic interaction usually favors a compression 
of the lamellar complex - at least if the underlying as- 
sumptions of our model are fulfilled. These assumptions 
are discussed in Sect. ^ where we also present some con- 
clusions. 



2 Free energy of model lipoplex 

The aim of the following calculation is to determine the 
electrostatic contribution to the undulations of a lamellar 
stack of membranes with DNA molecules intercalated in 
between. Our model system consists of two constituents, 
the membranes and the DNA molecules. The membranes 
have a uniform thickness t and carry positive charges on 
both sides. The surface charge density is given by a/2 
on each side of the bilayer and is assumed to be uni- 
form. In our model the membranes are perfectly trans- 
parent for the electric field lines, i.e., we have a homo- 
geneous dielectric constant throughout the lipoplex. The 
bilayers are flexible with a bending rigidity kc- The DNA 
molecules are modelled as infinitely long rigid rods of ra- 
dius r. For simplicity, we assume the negative charges of 
the DNA molecules to be located along their middle axis 
with the linear charge density —p. Following the experi- 
mental observation of a lamellar stack with DNA forming 
smectic arrays we arrange the components of our model 
in the following way (cf. Fig. 1). All membranes are par- 
allel to the Xy-plane with their midplanes at the posi- 
tions z = 0,±2{r + t/2) ,±4:{r + t/2) ... The DNA rods 
are aligned in the y-direction. The interhelical spacing 
between neighboring DNA molecules is constant and is de- 
noted hy d = 21: 1 q. Excluded volume requires that d > 2r. 
The rods in one layer are located at a; = 0, id, ±2d..., in 
the neighboring layers they are displaced by (i/2, i.e., they 
are at the positions x = ±d/2, ±3d/2, ... etc. Furthermore, 
the rods are assumed to be always attached to the two 
neighboring membranes. 
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unit cell 



Fig. 1. Schematic view of the model lipoplex for the case of 
flat membranes (see text for detail) 



The electrostatic interaction between the charges is 
calculated within the Debye-Hiickel approximation. In this 
approximation the potential ^ is determined by A(l> — k^(1> 
with the appropriate boundary conditions. Here k^^ de- 
notes the Debye screening length that is given hy k^^ — 
{SnUslB) where Ug is the bulk salt concentration and 
Ib — e^/ekBT is the Bjerrum length (e is the unit charge, 
fcsT is the thermal energy and e the dielectric constant; 
e w 80 in an aqueous solution). The total electrostatic 
contribution F^i to the free energy of the system is given 
by the sum of the (screened) electrostatic interactions and 
the translational entropy of the counterions 



1 



dSa'<P 



(1) 



The integration extends over all charged surfaces of the 
system with a' being the corresponding charge densities. 

We ask the following question: How are the membranes 
deformed by the electrostatic interaction? In order to an- 
swer this question we will calculate the induced undula- 
tions of the membrane up to the first order in the defor- 
mation amplitude. 

Consider the membrane at z = 0. Electrostatics in- 
duces a deformation z = u (x) around the flat state, z = 0. 
Due to the symmetry the deformation profile is of the form 



u {x) = ^^fln cos (nqx) 



(2) 



where the hat denotes summation over odd n only. This 
undulation leads to the following curvature energy fbend = 



(kj2) {V^uY (per area A) : 



kc sr-^ 2 4 4 

fbend = -^2^ a^n q 



(3) 



In order to calculate the electrical free energy, Eq. we 
compute first the electrical potential ^'^^^ {x, z) induced by 
the charges on the upper surface of the membrane (note 



that <P^^j^ is translational invariant in ^-direction) 



At that 



charged surface, i.e., at z = t/2+u (x), we have the bound- 
ary condition 8'^'"^^^ / dn = —l-Kaje which is here of the 



form 



y^,o-nnq sin {nqx) 



(") 

M 



dz 



-'-^ (4) 



(up to terms of the order a^J. By expanding up to 
first order in the amplitudes a„ we find the following form 
of the potential above the membrane (z > i/2 + u (.x)) 



ix,z)^ip("> (a;,z) + > W") ix,z) 



(") 



(5) 



Each (p*^"' fulfills the Debye-Hiickel equation separately. 
They can be expanded in Fourier series ip^"'^ = -Bm (z) 
cos (mnqx) where i?m (z) — bm^ exp (— Knmz) with = 

' -|- (nmq)^ . The coefficients 6m'' follow from the bound- 
ary condition at the membrane together with the fact that 
due to symmetry "P^^^ (x, z < t/2) = ^^"^ H'^/q) — x,t ~ z). 
We find that only the coefficients 5^"^ are non-vanishing 

and are given by fe^""* — 7rf7K/eK„. This leads to (for z > 
t/2): 



en L 



-K„(z-t/2) 



cos {nqx) 



(6) 



The total potential <1>m induced by the membrane at z = 
u {x) is the sum of the contributions of the upper charged 

boundary, cf. Eq. ||, and of the lower one, <Pm — 



(I) 



j^j. Using <1>^1\ (x, z) — <P)^/ (x, z + t) we find: 

27rcr ' 

en 



<Pm (x, z 



cosh 



Ea„ — cosh 



' COS (nqx) 



(7) 



Furthermore, the potential induced by the line charge of 
the rod has the form 



(i?) = -^Ko{kR) 



(8) 



where R is the distance from the line and Kq is a modi- 
fied Bessel function with Kq (x) ~ — In x for x <^ 1 and 
Kq (x) ~ (7r/2x)^^^ exp (— x) for x ^ 1. The total electri- 
cal potential ^ is the sum of the potential <Pi that follows 
from all membranes and the potential that is due to 
all the rods: ^ = + ^2- 

We calculate now the total electrostatic contribution 
to the free energy per unit cell. A unit cell has the width 
d (in AT-direction) and a height that corresponds to the 
(average) distance between neighboring layers (for the case 
of a flat membrane - depicted in Fig. 1 - this height equals 
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h>0 




h<0 




Fig. 2. Membrane undulations in lipoplexes. Shown are the 
two cases (a) ^ > 0: compression and (b) h < 0: swelling of the 
lamellar stack 



2r+t). According to Eq. |], we obtain the total electrostatic 
energy (per unit cell) by integrating the total potential 
over all charged surfaces that lie within this cell. The unit 
cell in Fig. 1 contains three charged surfaces, 5'^/\ S'^'f' 
and Sr. S"^"-* is a stripe of the upper surface of one bilayer 

that is uniformly charged with the density a/2. s'^J] is 
the corresponding lower charged surface. Sr is the surface 
carrying the charges of one rod; we assume this to be the 
surface of a cylinder with radius 6r <^ r and charge density 
—an = —p/27T5r. It follows that the electrostatic energy 
(per area) has three contributions: The inter-(and intra-) 
membrane interaction /m, the membrane-rod interaction 
Jmr and the interaction between the rods /r. Thus 



fel = /m + .fMR + fR = ^ I dS-Pi 



2A 



Sm 



dS<p2 



(9) 



Sr 



Here we made use of the identity — dSaR<Pi — dS 

We start by calculating the change of the membrane- 
membrane interactions induced by their undulations 
(up to first order in a„ ). The position of the midplane of 
the kth membrane is given by 



with fc = 0, ±1, ±2, ... and h = J^n'^n. Fig. 2 shows schematic 
views of lamellar structures that are compressed ~ case (a) 
with h > - and swollen - case (b) with h < 0. Denote 

(k) 

the contribution of fcth membrane to the potential by <P\^j . 
Then 



Tm 



qa_ 

in 



E 



fc=-c 



27r/g 



2 /Kt 



a„ cos (nqx) 



TTcr^ cosh (^) cosh (nr) ttct^ cosh^ (^) J2n'^n 
~ smh{K{r + t/2)) ^ esinh^ (^(r + I)) ^ ' 

Equation |ll| shows that a swelling of the system {h = 

X^n'*" ^ 0) decreases the membrane-membrane interac- 
tion whereas a compression (h > 0) is unfavorable. Note 
that we neglected terms of second order in the a„. As 
can be seen from Eq. |^ terms of the form lead to a 
renormalization of the bending constant, k'^ = kc + Skl^^. 
It can be shown that dk\^i = Sttct^/Se/c^ w T/k^IeXoh 
(Xgc = e/2-KlB(T is the Gouy-Chapman length) p^, ^ , |3^ , 
^,^|3^,^,^,|4^. For a wide range of parameters one 
has Sk\^i <^ kc- In the following we use the bare bending 
rigidity fcc, keeping in mind that it has to be replaced by 
k'^ when Sk\^i is comparable to kc. 

We estimate now the contribution of the membrane- 
rod attraction. Consider the rod at the position x = 
and z — —r — t/2 + h. The rod is located within an in- 
finite stack of membranes. This can be accounted for by 
simply summing twice over the contributions of all the 
membranes that are located above the rod, i.e. /mr = 
- C^P/d) E^o (0. -r-t/2 + h). From Eq. | follows 
that <p']!^j (0, z) is given by (note that z < 0) 



27rCT / Kt 

cosh [ — \ e 

en \ 2 



K,{z-2k(r+t/2-h)) 



^ ^cosh(^)e'="(--2K-+l)) (12) 



(-1) 



fe+i 



The contribution of the first term of Eq. 12 to Fmr is of 
the form: 



' MR. 



apq cosh(Kt/2) 
en sinh (k (r -|- t/2)) 



crpg coth (k (r I)) , f Kt\ 

— 7-4 ^ cosh — > a„ (13) 

e sinh(K(r + I)) V 2 / ^» 

The second term of Eq. |l^ leads to the following expres- 
sion: 



(2) _ apqn 

J MR — 



X ^ cosh(K„i/2) 

cosh («:„ (r -|- 1/2)) " 



(14) 



The total membrane-DNA contribution /m_r = Z^^'' 



(2\ 

J MR fervors a compression of the lamellar stack - thus 
- — - constituting a competing mechanism to the mcmbrane- 

1) > a„ cos (nqx) (10) membrane repulsion. 



MR 



Uk (x) =2k{r + t/2~h) + {- 
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We are left with the calculation of the interaction en- 
ergy between the rods. We focus here one two important 
cases. Case 1: Kt ^ 1 and nd <^ 1 ("vertical screening"): 
In this case the interaction between rods in different layers 
is negligible compared to the rod-rod interaction within 
the same layer. Then it is sufficient to sum over the con- 
tributions of all rods to the left and to the right of the 
given rod : 



2 °° 2 

JR = — Ko [ndk) 



ne 



k=l 



ire 



dkKo (ndk) 



Attek 



(15) 

Case 2: Kt <^ 1 and Kci <C 1 (weak screening): In this 
case all rods contribute to the interaction energy. After 
some algebra we arrive at 



9 9 

p q 



2 2 

p r 



ATTEK^ir + t/i) 47r£«;2(r + i/2)2^n " 



(16) 



As expected, in both cases the repulsive rod-rod interac- 
tion favors swelling. 



3 Undulations in isoelectric complexes 

We consider first the lamellar complex in equilibrium with 
a solution of free DNA strands. We ask: What is the inter- 
axial distance d between the DNA strands in the lipoplex 
that minimizes the electrostatic free energy of the com- 
plex? As pointed out by Bruinsma [p^ counterion release 
will control d; here, however, we determine the equilibrium 
spacing for the case when there is no counterion release, 
i.e., we assume the rods being below the Manning thresh- 
old. We also neglect entropic changes due to the adsorp- 
tion of free DNA strands into the lipoplex. We will show 
that, as a result of the geometry, such a system will equi- 
librate at the isoelectric point - if the electrostatic inter- 
action is sufficiently long-ranged. In the following we only 
account for the contributions independent of the a„'s and 
treat the contribution of membrane bending afterwards as 
a perturbation. 

Let us first consider the case of high ionic strength 
where nr ^ 1 (strong screening). Then the free energy 
per area is given by 



Anape 
end 



(17) 



i.e., by the membrane-rod attraction, Eq. |l^; other terms 
are negligible. It follows that the minimum is at d 0. 
Excluded volume interaction between the rods will lead 
to d = 2r. Clearly, in the case of strong screening as a 
result of the short range of the electrostatic interaction 
the lipoplex is equilibrated far from the isoelectric point. 
The resulting complex is "overcharged" by the DNA rods. 
(A similar situation occurs for the adsorption of rods on 
an oppositely charged surface, cf. Ref. Q). 

We discuss next the two cases introduced above. Case 
1: Kt ^ 1 and Kd <^ 1 (vertical screening): From Eqs. ^ 



and we find (up to terms of the order Kd) 



fu 



2lTap TT/O 



end end^ 



(18) 



ftot is minimized for d — diso — p/a which corresponds 
to the isoelectric point of the complex, i.e., the point at 
which the charges of the cationic lipids and of the DNA 
are exactly balanced. Case 2: nd 1 and nd <^ 1 (weak 
screening): From Eqs. O and n6^ follows 



2'K(7p 



np 



eK'^{r + t/2)d eK^{r + t/2)d^ 



(19) 



Again the free energy is minimized at the isoelectric in- 
terhelical spacing d = diso — p/c- 

Thus in the limiting case k ^ (no salt, no screening) 
the lipoplex is forced to be at the isoelectric point. This 
is an artefact of our model which does not account for 
counterions that would be present in lipoplexes with an 
excess of DNA molecules or cationic lipids. This holds for 
both cases, for Case 1 corresponding effectively to decou- 
pled two-dimensional layers and for Case 2 which is truly 
three-dimensional. Therefore our theory is only applicable 
to isoelectric lipoplexes where all counterions are expected 
to be released. 

We consider now the undulations occuring in lipoplexes 
in general and then focus again on the isoelectric point. 
The change of the total electrostatic free energy as a func- 
tion of the deformation follows from the Eqs. 



;yas a tunc 



4^^^ ) «2(r+|)' 



for > 1, < 1 

for < 1, Kd < 1 
(20) 

The change of the total free energy due to bending is given 
by the sum of the electrical contribution Z\/ei and the 
bending energy jhend-, Eq. ^ Minimizing ]tot with respect 
to the amplitudes a„ leads to a„ ~ Ajn^ with 



A 



for Kt > 1, Kd < 1 
for Kt < 1, Kd < 1 



(21) 



Thus the deformation modes decrease rapidly with in- 
creasing n. Now we are in the position to calculate the 
deformation of the membranes. Inserting a„ = A/n'* into 
Eq. ^ we find the following deformation profile 

cos (nqx) 



u(x) ^ A 



TTX 

~V1 



/ '^r-^ = ^ \ 

'^n n4 12 12 48 

(22) 

where the polynomial expression is valid for < < tt 
(the continuation outside this interval follows from the 
symmetry of the configuration). A good approximation 
for all values of x (relative error smaller than 1.5%) is 
given by u [x) — A cos (ga:). 

The coefficient A at the isoelectric point of the complex 
21X0 1 q = p is given by 



ekcq' 





:r for K^ 3> 1 , Kd ^ 1 
for Ki < 1, Kd <C 1 



(23) 
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In the case of vertical screening we find a positive (and 
K-independent) value of A and thus a positive value of h, 
h — tt'^A/90, corresponding to a compression of the iso- 
electric lamellar stack. Interestingly, in the case of weak 
screening the undulations disappear. In fact, as long as the 
vertical screening is operative the membrane-membrane 
repulsion is smaller than the membrane-rod attraction in 
the isoelectric lipoplex and therefore we find a compres- 
sion of the lamellar stack. For weak screening, the rod-rod 
repulsion between different layers cancels this net attrac- 
tion, cf. Eq. Let us consider typical values for p, cr, 
e and kc, say p = ejX.lk (DNA), a = e/lOOA^, e = 80 
(water) and fee = 20kBT. For these values we find (in the 
case of vertical screening) h w lA, i.e., the undulations 
are rather small. 

Finally, we estimate how the undulations of the iso- 
electric complex disturb the interaxial spacing between 
the rods and in turn move the complex away from its iso- 
electric point. We consider the case of vertical screening 
(Case 1, Kt ^ 1, Kci <C 1). In this case the amplitudes of 
the undulations depend strongly on the interhelical dis- 
tance, namely A ~ d^, cf. Eq. Inserting h — X^n'^" ~ 
7r''A/90 into Eq. EG (k^ :$> 1, nd 1) we find two correc- 
tion terms to Eq.llq, namely n'^a'^d'^/ (l440£^fcc) from the 
membrane-membrane repulsion and —n^a^ pd^ / (720e^fcc) 
from the membrane-rod attraction. Evidently, the mem- 
brane-membrane repulsion favors smaller values of d that 
lead to smaller undulations, whereas the membrane-rod 
interaction is enhanced for larger undulations, i.e., larger 
values of d are favorable. At the isoelectric point d = p/a 
the correction term from the membrane-rod interaction 
exceeds the other term and as a result the interhelical dis- 



tance is slightly increased, d 



Ad- 



Ad, with 



Kp 



720 ekra^ 



(24) 



In that sense the undulations lead to an effective repulsion 
between the DNA strands proportional to d^ (as long as 
one is close enough to the isoelectric point). Undulations 
are one of several mechanisms that might be responsible 
for the increase of the interhelical spacing with increasing 
salt concentration as it is observed experimentally (cf. Ref. 
H for details). 



4 Discussion 

The main idea of the preceding analysis is to give a simple 
estimate of the role of the electrostatic interactions in a 
lamellar lipoplex. In many instances some of the underly- 
ing assumptions are not fulfilled. But even in this case our 
model might give an idea about what the contributions of 
the electrostatics to the overall conformation might be. 

One severe approximation is the assumption of a trans- 
parent membrane. The lipid bilayer represents a low di- 
electric slab that - depending on its thickness - might 
screen most of the electrical field so that charges (say of 
phosphate group on a DNA molecule) are not "seen" on 
the other side of the membrane. As a rule of thumb, a 



membrane that is much thinner than the distance D of 
a charge from the membrane might be considered to be 
transparent for this charge whereas a thicker membrane 
(of thickness t > D) is opaque and can be approximated 
by an infinitely thick slab. In that case the effect of the 
low dielectric lipid can be accounted for by the use of an 
appropriate image charge - which is a simple task for a flat 
membrane but difficult to handle for an undulating one. 
The two cases t < D and t > D are elaborated in some 
detail in Footnote 2 in Ref. The typical thickness of 
a lipid membrane is 24A which is of the order of the di- 
ameter of the DNA rod (20 A), i.e., one is in the crossover 
regime between the two cases. In any case, the presence 
of the low-dielectric lipid will lead to a modification of 
the simple situation discussed in this paper. It should be 
expected that the partial confinement of electrical field 
lines emanating from the DNA rods by the neighboring 
membranes favors a swelling of the lamellar stack. 

Another feature not considered in this study is the 
demixing of neutral and charged lipids within the bilayers. 
There might be at least three effects, (i) The electrostatic 
attraction drives the cationic lipids towards the DNA rods, 
resulting in a depletion of charged lipids in the mem- 
brane parts in between two neighboring rods in the same 
layer. This reduces the membrane-membrane repulsion be- 
tween neighboring membranes resulting in an increase of 
the compression of the lipoplex. This effect should be im- 
portant if the average mole fraction of cationic lipids in 
the bilayers is low. (ii) For the opposite limit of highly 
charged membranes their surface charge density a might 
exceed the surface charge density p/2'Kr of the rods. In this 
case an enhancement of neutral lipids close to the DNA 
is expected that allows for a better matching of the two 
charge densities. These two effects were indeed observed 
in the numerical study by Harries et al. |^ (cf. Fig. 7 
in that paper). [Hi) Finally, membrane undulations may 
also affect the charge densities on each side of the bilayer 
and vice versa. A depletion of cationic lipids on one side 
of the membrane will be accompanied by an enhancement 
on the other side due to the symmetry of the arrangement 
of DNA rods, cf. Fig. 1. This might lead to a spontaneous 
curvature of the bilayer of either sign which in turn af- 
fects the membrane undulations discussed in our study. 
It is clear that the competition of these three effects can 
lead to a rather complex behavior of the charge density 
profile of the lipids along the AT-direction. In our model 
we do not account for these effects. We expect that our 
approach has to be modified especially when a large frac- 
tion of the lipids is neutral. In that case the description 
of local, highly charged pinches |^ might be more appro- 
priate (cf. the discussion of the polypeptide lipoplex [ pO) 
given in the Introduction of our paper). 

Recently it has been possible to non-ambiguously de- 
tect undulations in a lipoplex|Q. From a careful analysis 
of the data on the lipoplex presented in Ref. if was 
possible to construct an electron density map revealing 
its high resolution structure. The structure shows undu- 
lations with an amplitude of a few Angstrom leading to a 
compression of the lamellar structure as depicted in Fig. 
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2a. Furthermore, the amplitudes show a sharp increase for 
larger inter hehcal spacings. Both observations arc in qual- 
itative and semi-quantitative agreement with the resuhs 
of our model calculation, cf. Eq. It is worth noting that 
the lipids of these lipoplexes form a lipid-gel phase [|| . In 
this phase the lipid bilayers show a high compressibility 
allowing the DNA-induced deformations to cross nearly 
unperturbed through the bilayer (as implicitely assumed 
in our model) . The induced undulations might be the pre- 
vailing mechanism for the interlayer coupling that leads to 
the rectangular columnar superlattice of the DNA strands 
observed for this class of lipoplexes. 
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